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I. INTRODUCTION 

Dispersive shock waves (DSWs) are observed in non- 
linear optics in systems described by the nonlinear 
Schdrdinger equation, when the so-called hydrodynam- 
ical reduction is valid [THE]. The introduction of a small 
amount of disorder competes with nonlincarity and ham- 
pers the shock formation [8] . This makes the DSWs 
an appealing framework to study the interplay between 
randomness and nonlinear waves, a subject of growing 
interest as witnessed by recent theoretical [OHTTj and ex- 
perimental studies [TMT4] . 

At variance with the ordered systems [15], the direct 
observation and characterization of optical shock waves 
in the presence of structural randomness is burdened by 
several technical difficulties in identifying an appropri- 
ate nonlinear medium, feasible excitation conditions, and 
relevant observables. There are several possibilities to 
characterize the excitation of undular bores and related 
phenomena |16H21| : the very definition and observation 
of the wave-breaking phenomena in the presence of dis- 
order and nonlinearity is an open issue. This calls for 
an extensive development of experimental techniques and 
the use of multiple methods to characterize the DSWs. 

In this paper we give a detailed review of our exper- 
imental investigations of the hydrodynamical regime in 
the generation of optical shocks during nonlinear opti- 
cal propagation in a thermal defocusing nonlinearity of 
a continuous wave (CW) laser beam. The hydrodynami- 
cal regime is achieved when the nonlinear length is much 
smaller than the diffraction and losses (absorption and 
scattering) lengths. 

Our experimental technique allows the direct obser- 
vation of a propagating initially Gaussian laser beam 
in a thermal nonlinear liquid with controllable disorder 
obtained by a colloidal dispersion with a low index of 
refraction contrast. We show that, by increasing the 
strength of the nonlinearity, shock formation is enhanced, 
while, on the other hand, the random scatterers limit and 
ultimately inhibit the wave-breaking phenomenon. We 
quantify such a competition by analyzing both the laser 
beam along the propagation direction and its far field dis- 
tribution intensity: these allow the measurement of the 
relevant scaling laws that relate the shock position [7] [55] 



and the post-shock wave-vector spectrum with the input 
beam size, power, and the scale and strength of disorder 
0. 

These observables, namely the shock formation point 
and the output wave-vector spectrum, exhibit a thresh- 
old and determine a phase diagram identifying parameter 
regions where the shock occurs. 

The paper is organized as follows: first we briefly re- 
view the theoretical framework. We then illustrate the 
experimental setup and the characterization of the sam- 
ples used in the experiments. We then review the re- 
sults in three different sections, corresponding to the ex- 
perimental characterization of different observables. We 
present the analysis of the beam intensity distribution 
both during propagation and exiting the samples. We 
conclude with a summary section of the obtained results. 

II. THEORETICAL FRAMEWORK 

In our experiments, CW laser beams propagate in dye- 
doped dispersions of dielectric colloidal beads. The beam 
is partially absorbed and scattered, activating the inter- 
play between thermal-defocusing and spatial-disorder. 

Neglecting, in a first approximation, the spatial nonlo- 
cality [23] . the refractive index perturbation in the pres- 
ence of nonlinearity and disorder to the bulk index uq is 
written as: 

An = n 2 I + An R (X,Y,Z) (1) 

where n 2 < takes into account the considered defo- 
cusing Kerr effect, I is the optical intensity and An^ 
represents the random perturbation due to the colloidal 
beads. 

The propagation of a TEMoo Gaussian beam inside 
the medium is described by the paraxial wave equation 
for the complex envelope, A, of a monochromatic electric 
field E = (^^Y^AexpiikZ - iuT), 

2ik £ ^ A +V 2 XY A + 2k 2 —A = 0, (2) 
oZ no 

where k = 27rno/A is the wave- vector, c the velocity of 
light, and eo is the electric permittivity of free space. In 
Eq. (pi) A is normalized such that / ~ |A| 2 . Indicating 
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with Iq the input peak intensity, wq the input beam waist, 
L n i = no/ (feo 1^2 1 Jo) the nonlinear length scale, and intro- 
ducing the scaled coordinates x, y, z = X/wq, Y/wq, Z/L, 
and the normalized field ip = A/y/Io, we obtain the fol- 
lowing dimensionless equation: 

fc^ + fv^-MV + W = 0, (3) 

where e = L n i/L = ^L nl /L d , being L = \/L n iL d and 
Ld = kwQ the diffraction length, and Ur = Aur/ (n2lo)- 
The quantity e measures the strength of the nonlinear- 
ity with respect to the diffraction: a small value for e 
implies negligible diffraction and a pronounced nonlinear 
response. Ur is the ratio between the perturbation of 
index due to the disorder and the nonlinearity. 

Setting in Eq.|3| ip = yj p(r, z) exp[i^>(r, z)/e] and re- 
taining only the leading order in e, we obtain the follow- 
ing equation for the phase <f>: 

<t>z + \{<t>l + <t>l) + P-U R = Q (4) 

Limiting to one dimension (ID, d y — 0), performing 
the transverse derivative of Eq.Q, and defining a velocity 
field equal to the phase chirp, u = 4> Xl we have 

u z + uu x + d x (p- Ur) = 0. (5) 

In the homogeneous case ip =const.) and for an ordered 
medium (Ur = 0), Eq. ^ takes the form of the Hopf 
equation [5] , the solution of which can develop disconti- 
nuities in the velocity profile, u x — > oo, and hence gives 
rise to shock waves. 

Here we remark that from the hydrodynamical approx- 
imation, a threshold in the nonlinearity is present: in 
fact, the approximation holds true when L n i -C Ld- An- 
other threshold arises from the term Ur — Anfl/(n 2 ^o)j 
corresponding to the existence of a critical value for the 
amount of randomness, above which it is expected that 
no shock occurs: when the random index perturbation 
Ahr becomes comparable with the nonlinearity ^Iq, the 
material refractive index fluctuations are so pronounced 
that the nonlinear effect is totally masked. 

Correspondingly, in our experiments (see below) in ab- 
sence of disorder, we find a threshold in the laser power, 
while in the presence of disorder a threshold emerges also 
in the amount of randomness. 



III. EXPERIMENTAL SETUP. 

In the hydrodynamic limit, DSWs are expected to 
occur when the nonlinearity is dominant compared to 
diffraction; nevertheless the diffraction, which is initially 
negligible, starts to play a major role in the proximity of 
the wave-breaking point, and regularizes the singularity 
by means of the appearance of characteristic oscillations 
(undular bores). 



Besides these regularizing oscillations, the singularity 
in the field phase and amplitude also results in a diffrac- 
tion enhancement, evident in the funnel shape along the 
propagation direction (see below) appearing with the in- 
crease of the input power. This shows that the shock 
involves the spatial spectrum of the beam as detected in 
far field measurements. 

Near-field configuration - In the near-field configura- 
tion our setup [Fig. [ija)] allows a direct visualization of 
the propagating beam profile, i.e., the intensity as func- 
tion of the transverse coordinate, X, and of the prop- 
agation direction Z. This enables the identification of 
the shock point Z s as the propagation distance at which 
the maximum chirp occurs (see below) . Typically, a CW 
laser at wavelength A = 532nm is focused inside the sam- 
ple. The beam waist in the focus is wq = 10/im. The 
near-field configuration is sketched in Fig. [ija). A 1cm 
x 1cm x 3cm glass cell is used and the laser beam prop- 
agates along the 1 cm side. Top images of the fluores- 
cence emission are collected by a MZ16 Leica microscope 
placed perpendicularly to the propagation direction, Z, 
and recorded by a 1024 x 1392 pixels CCD camera. 

Far-field configuration - In Fig. [ljb) we show the setup 
for the far-field measurements. The CW laser beam is 
focused inside the sample (wq = 50/im) . The liquid sam- 
ples are placed in a 1mm xlcm x3cm glass cell, the laser 
beam propagates along the 1mm side, the cell is placed 
in vertical direction in order to moderate the effect of 
heat convection. As shown in Fig. [ijc) the intensity 
distribution of the Fourier transform of the transmitted 
beam is collected by a CCD camera placed at the focal 
length from the collecting lens. We calibrate the CCD 
detector by fitting with the Airy function the experimen- 
tally obtained Fourier transform of a 500/im diameter 
pinhole, placed on the exit face of the cell. The angu- 
lar spreading is related to the transverse wavevector as 
(27r/A)sin(0). 

Figure [TJd) shows the mutual positions of the focus 
plane, the shock plane and the output plane. 



IV. SAMPLE CHARACTERIZATION 

As in previous experimental works, we use the thermal 
Kerr-like defocusing nonlinearity of absorbing dye-doped 
liquid media [3 H • 

Our samples are aqueous solutions of Rhodamine B 
(RhB). We tailor the degree of absorption and nonlinear- 
ity by varying the concentration of RhB (cr^b) from 0.05 
to 0.2mM. We add disorder by using monodisperse 1/zm 
diameter silica (SiC^) spheres. The degree of random- 
ness is fixed by varying the concentration of Si02 (csi0 2 ) 
from 0.005 to 0.04w/w, in units of weight of silica parti- 
cles over suspension weight. In terms of refractive index 
perturbation, the amount of disorder can be estimated 
by the following relation: 

(A"/j) 1/2 = c Si o 2 PH 2 o(nsi02 ~ nH 2 o)/psi0 2 , (6) 
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FIG. 1: (Color online) Experimental setup: (a) detection of the top fluorescence emission of the beam; (b) configuration for 
the collection of the far field intensity; (c) details of the optical setup of panel (b); (d) sketch showing the position of the shock 
plane inside of the sample. 



being n Sl o 2 (nfJ 2 o) and p S io 2 (ph 2 o) the refractive in- 
dex and the density of the Si0 2 (H 2 0), respectively. 
The angular brackets in Eq. ^ denotes volume aver- 
age. Being the silica (water) density psio 2 = 2g/cm 3 
{ph 2 o ~ lg/cm 3 at 25°C), for the considered range of 
c$io 2 concentration, (An^) 1 / 2 varies between 4 x 10~ 4 
and 32 x 1CT 4 . Therefore, since from the theory a thresh- 
old in the disorder amount is predicted when (An 2 }) 1 / 2 
becomes comparable with the nonlinear perturbation 
l n 2|A) — 10~ 3 , such a threshold is expected for the silica 
concentration csi0 2 — 0.030w/w as it was confirmed by 
our experiments (see below). 

In our samples there are two leading loss mechanisms: 
(ii) absorption due to the RhB dye, and (ii) scattering 



due to SiC>2 particles. We find that scattering losses are 
predominant; this is shown in Fig. [2] where we compare 
the images of the transverse beam intensity distribution 
versus the propagation direction Z at two different SiC>2 
concentrations, fixed laser power P and RhB concentra- 
tion c RhB , 

Figure [2ja) shows the top fluorescence of the laser 
beam in a pure dye sample (csio 2 = Q)j Fig- pHb) gives 
the corresponding far-field. Figures 12] (c) and u\d) re- 
port the case of a silica-dye sample at csio 2 — 0.03w/w; 
the beam is more diffused and the far field reveals an 
enhanced spectral content. In the presence of disorder 
( c Si0 2 > 0) the transverse spread of the beam along Z 
is enhanced. This is clarified in the analysis reported in 
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FIG. 2: (Color online) (a)-(c) Near-field images of the fluores- 
cence emission of the propagating beam at fixed laser power 
(P = 8mW) and two different particle concentrations, (a) 
csi0 2 = 0, (c) csi0 2 — 0.03w/w; (b)-(d) corresponding inten- 
sity profiles of the far-field images of the transmitted field at 
the exit face of 1mm thick cell. 
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FIG. 3: (Color online) (a)-(b) Intensity profiles taken at 
three different Z position in the images of Figs. [2^a)-[2jc): 
dashed line corresponds to Z — 0.3mm, dot-dashed line to 
Z — 0.7mm, and the continuous line to Z — 1.0mm. (c) Ex- 
ponential (quasi-linear) decays of beam intensity calculated 
from images of Fig. |2fa) (dashed line) and Fig. ^c) (dot- 
dashed line). 



Fig. ^a) and[3]Jb), which show the intensity profiles at 
three different Z positions for the O.lmM pure dye solu- 
tion and dye solution with silica at csio 2 = 0.03w/w con- 
centration, respectively. At variance with linear absorp- 
tion, scattering due to SiC>2 beads, broadens the beam. 
Figure |3jc) shows the average intensity Vs Z calculated 
from the images of Figs. [2]ja) and [2jc) and the expo- 
nential decays that fit the data. The fitting coefficients 
of the exponential decays give: the absorption length, 
Labs = 1.6mm, and the losses (i.e., absorption and scat- 



tering) length, Lies — 1.2mm, for the pure dye solution 
and for the 0.03w/w silica-dye solution, respectively; this 
implies that the effect of the particles on the losses is very 
small. 

In summary the above analysis shows that the role 
of disorder is predominantly to introduce random phase 
modulation. Moreover, since absorption does not qualita- 
tively affect shock formation, the disorder induced phase 
scrambling is predominant over all the loss mechanisms 
in determining the shock point Z s measured below. 

V. SHOCK POINT 

In this section we report the procedure to identify the 
shock point Z s and to determine the threshold for the 
wave breaking in terms of laser power and Si02 concen- 
tration. 

In Figs. gaj-gc) andga) and^c), we show the im- 
ages of the propagating beam versus Z direction at low 
and high laser power, respectively. At low power, i.e., 
P = lOmW, no nonlinear effects are visible. In Figs. 
|4^d)-|4jf) and |4jd)-|4jf) we show the corresponding im- 
ages of the output intensity field. 
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FIG. 4: (Color online) Top panels: low power (P = lOmW) 
images of the fluorescence emission of the beam along Z at 
different concentrations of silica spheres: (a) csi0 2 — 0, (b) 
Csi0 2 = 0.017 w/w, (c) c$io 2 = 0.030 w/w; bottom panels: 
corresponding images of the transmitted (in the X-Y plane) 
intensity at the exit facet of 1mm cell. 

Conversely at higher beam power both the effects of 
nonlinearity and of disorder are evident simultaneously. 
The two effects are competing as evident by the shock 
features, i.e., an augmented beam diffraction and the ap- 
pearance of the undular bores, are enhanced by the laser 
power and inhibited by the Si02 concentration, clear 
from the transverse and longitudinal intensity profiles of 
Fig. [5} 

We determine the shock point Z s from the intensity 
profile of Fig. [4] and [5j recalling that the shock is origi- 
nated from a singularity in the phase chirp \d(f>/dX\ — > oo 

01221 E3I- 

To retrieve the phase singularity from the intensity 
profile we used the following argument: in the hydrody- 
namical approximation the laser beam is mainly affected 
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FIG. 6: (Color online) scattered dots are the calculated steep- 
ness along the Z direction for three different powers. The solid 
lines are polynomial fit of the steepness curves to identify their 
maximum value, indicated by the arrows. 

by the defocusing nonlinearity, in a regime of negligi- 
ble losses and diffraction. Hence at first approximation 
the phase is proportional to the refractive index pertur- 
bation, which in turn depends on the intensity profile 
because of the Kerr nonlinearity: 

4>(X,Y,Z) = ^An[I(X,Y,Z)}. (7) 
n 

From Eq. Q we can estimate the occurrence of the sin- 
gularity in the phase from the intensity profiles, in fact: 

Vx,Y(b{X, Y, Z) oc Vx,yI{X, Y, Z). (8) 

Equation ffify shows that the point of maximum phase 
chirp is given by the maximum derivative in the inten- 
sity profile, this allows the estimation of the shock point 
as follows: we calculate the transverse derivative of the 
intensity normalized to the peak value, 7jv, and we define 
the steepness S(Z) as the maximum with respect to the 
transverse coordinates of such a derivative: 
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FIG. 7: (Color online) (a) measured Z a vs P for different 
csi0 2 i (b) power-disorder diagram from the propagation mea- 
surements: filled circles are the threshold power calculated 
from (a), dashed line is a boundary due to the experimental 
available observation window, and the dot-dashed line is the 
boundary as estimated by the theory. 

The shock point, Z s , is finally defined as the position of 
the maximum steepness versus Z: 

Z s = max z [S(Z)}. (10) 

In Fig. [6] we show the steepness curves S(Z) at three 
different laser power P. The point of shock occurs at 
propagation distances consistently smaller with the in- 
crease of the incident laser beam power, a signature that, 
for a fixed level of disorder, the increase of the nonlinear- 
ity enhances the shock formation. We note that the curve 
corresponding to the lowest P shows a monotonous trend 
and reaches its maximum value at the edge of the obser- 
vation window. This implies the existence of a thresh- 
old value of P below which Z s assumes a constant value 
(equal to the size of the observation window Lq ~ 1mm) . 

In Fig. Qa) we plot the calculated Z s vs P for all 
the prepared csio 2 concentrations. We observe that the 
threshold power at which Z s starts to decrease with re- 
spect to Lq becomes larger when increasing csio 2 > result- 
ing in a shift of the power threshold towards higher val- 
ues. In Fig. ^b) we map the threshold P in a disorder- 
power shock phase diagram. 

We remark that the obtained Z s values are in all the 
investigated cases always smaller of the absorption length 
L a b s , confirming that the absorption only marginally af- 
fects the shock formation that is instead connected to the 
phase scrambling due to the silica particles scattering. 



VI. INTENSITY CORRELATION AT THE 
SHOCK 

In the previous section we have quantitatively analyzed 
the top fluorescence near- field images of the propagating 
beam by calculating the shock position Z s . In order to 
also analyze the transmitted profiles (bottom panels of 
Figs. [5] and [6} we calculate the correlation function as 
follows: 



S(Z) = max XiY [V ' XtY I N {X,Y,Z)\. (9) 



C{P) = Ei d I P (iJ)Io(i,j)/EijIo(i,j)Io(iJ), (11) 
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FIG. 8: (Color online) (a) Correlation curves calculated from 
the far-field measurements calculated for the ordered sample 
and for three different increasing values of csi0 2 \ (b) power- 
disorder phase diagram as calculated from the transverse in- 
tensity correlation curves of (a). 



where P is the laser power, Ip(i,j) the intensity distri- 
bution on the CCD camera, with i and j pixel indexes, 
corresponding to a certain power P, and Io(i,j) is the 
reference image of the intensity distribution transmitted 
from the pure dye sample (cnhB = O.lmM) at the laser 
power P = 160mW, such reference image was selected 
as the first image clearly showing the post-shock rings. 
The function C(P) provides an estimation of the degree 
of coherence after propagation in the scattering samples. 

Figure [8]ja) shows the correlation curves C{P) calcu- 
lated from the bottom images of Figs. [5] and |6j we ob- 
serve that the curves grow up to a maximum value and 
then they start to decrease with the increasing power P. 
The power P of the peak value increases with the Si02 
concentration, meaning that in the presence of disorder 
a stronger nonlinearity is necessary to overcome the de- 
phasing effect due to the scattering with the silica parti- 
cles. Figure |8jb) shows the disorder-power shock phase 
diagram as calculated from the curves of Fig. |8ja) : the 
filled circles represent the threshold power P, defined as 
the power at which the maximum correlation between 
the ordered an the disordered samples is achieved. We 
stress that the shock phase diagram of Fig. [8jb) is in 
good agreement with that of Fig. [7](b) calculated from 
the Z S (P) curves; the slight discrepancy between the two 
phase diagrams derives from the different definition of the 
shock point. 



VII. SHOCK THRESHOLD FROM ANGULAR 
SPREADING MEASUREMENTS 

The characteristic post-shock annular structure and 
the diffraction enhancement displayed by the near-field 
transverse and longitudinal intensity distribution, reveal 
a non-trivial involvement of the wave- vector spectrum in 
the shock phenomenon. In this section we report the 
investigation on the far-field intensity distribution of the 
transmitted beam after 1mm propagation distance. Such 
an investigation allows us to measure the angular aper- 



ture 9. Fig. [9] provides a qualitative overview of the 
whole set of the far-field measurements. The panels on 
the left side [Figs. [9]^a)-[9]^c)] report images relative to the 
ordered samples {csio-, = 0) at fixed power P = 130mW 
and at various dye concentrations cr^b ranging from 0.05 
to 0.2mM; the right panels [Figs, ^a')-^')] re fer to the 
disordered samples at fixed power P = 140mW prepared 
at CRhB = O.lmM and varying csi0 2 between 0.005w/w 
and 0.038w/w. The way the nonlinearity and the disor- 
der affect the shock phenomenology, i.e., the appearance 
of the characteristic rings and the enlargement of the 
spectral content, reveals that their effect on the shock 
formation is opposite: the images [Figs^-Mc)] show 
an enhancement with the increase of crhb (i-e- f of the 
strength of nonlinearity) ; conversely those in Figs. [9ja,')— 
[9^c') show the inhibition of shock with csio 2 ■ 

Note that in images on the right the circular symme- 
try of the DSWs is lost because of the refractive index 
inhomogeneities. In other words, the shock wave has a 
partially randomized spatial distribution. In order to 
quantitatively analyze both the sets of measurements, 
we perform a radial average of the two-dimensional col- 
lected profiles and we estimate the angular aperture 9 as 
the full width half maximum since the profile appears as 
a single peak; and as the distance between the two lead- 
ing peaks when the profiles start to split because of the 
wave breaking due to the defocusing nonlinearity. 

In what follows we detail the results obtained for the 
ordered and disordered case [8]. 



A. Ordered case 

We study the occurrence of DSW in the pure dye so- 
lutions (csio 2 — 0) when varying input laser power P for 
different dye concentrations crkb- 

Figures [Ma)-[I(Jf ) display the collected images of the 
far-field intensity distribution when a low [Figs. |10[ a)- 
TO^b)] and high [Figs. [l0^d)-[To](f )] power laser beam im- 
pinges on the pure dye solutions. We note that the higher 
the dye concentration, the larger the spatial spectral con- 
tent due to the higher nonlinearity. 

Figure [Toj^g) shows the curves of the calculated angu- 
lar aperture 9 vs P as obtained for the different crkb 
concentrations. In these measurements both the control 
parameters contribute to strengthen the nonlinearity of 
the system. Consistently we find that, above a critical 
power, 9 starts to increase with P because of the speedup 
of the shock formation due to the augmented nonlinear- 
ity; the slope of the curves increases with c^Bi providing 
larger spectra at the same laser power P. 

Analogously to our analysis of the shock position Z s , 
we seek also for the angular aperture a threshold value for 
the laser power. Such a threshold power can be retrieved 
in the above mentioned power value, beyond which 9 
starts to linearly grow with P. We plot the threshold 
power values in the diagram of P versus cr^b m Fig. 

Mil). 
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FIG. 9: (Color online) Images of the far-field intensity distribution of the transmitted beam after 1mm propagation distance. 
The left panels group refers to ordered samples (i.e. fixed csio 2 — 0) at fixed power P = 140mW at various dye concentrations 
(a) crhb = 0.05mM, (b) O.lmM (b), (c) and 0.2mM. The right panels group shows the far field intensity profile in disordered 
samples (i.e., fixed cr^b = O.lmM) fixed laser power P = 130mW and at various SiC>2 concentrations (a') csio 2 = 0.005w/w, 
(b') 0.017w/w and (c') 0.038w/w. 
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FIG. 10: (Color online) (a)-(f): spectral content of the 
transmitted beam of the ordered samples (csio 2 = 0) at 
two laser power P and different dye concentration: (a),(d) 
CRhB = 0.05mM; (b),(e) O.lmM and, (c), (f) 0.2mM. Bottom 
panels: 6 vs P curves (g); threshold power P as calculated 
from curves of (g) vs Cijhs(h). 



B. Disordered case 

We consider the interplay between disorder and non- 
linearity in the DSW formation by dispersing the SiC>2 
particles in pure dye solutions at cr^b = 0.05mM and 
c RhB = O.lmM. Figures [TTJa)-(f), show the spectral 
profiles for different csio 2 an d laser power P at fixed 
CRhB = 0.05mM. At this dye concentration and any laser 
power P, no shock formation emerges from the spectra 
as can be retrieved also in the trend of 9 vs P reported 
in Fig. 11 (g). This is a signature of the fact that at 



the lowest prepared dye concentration the nonlinearity is 
counteracted by the disorder which prevents the appear- 
ance of any shock phenomenology. 
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FIG. 11: (Color online) (a)-(f): low dye concentration 
(criiB = 0.05mM) spectral content of the transmitted beam 
by disordered samples at two different laser power P and dif- 
ferent Si02 concentration: (a), (d) csi0 2 = 0.007w/w; (b), 
(e) 0.018w/w and (c), (f) 0.038w/w; (g) the angular aperture 
9 vs P: the curves show no threshold behavior in the presence 
of disorder. 

In Fig. [12] we show the case of the disordered samples 
obtained by the O.lmM pure dye solution. At the higher 
power P = 140mW the shock characteristic rings are 
clearly visible in the spectra corresponding to the lower 
concentrations of SiC>2, csio 2 — 0.007 and 0.018w/w and 
disappear at the highest concentration c$io 2 = 0.038. 

In Fig. 12 'g) we show the curves vs P calculated from 
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the images of the upper panels [Figs. [I2[a){l2[f)]. We 
retrieve the expected threshold behavior with respect to 
both the control parameters P and c$io 2 > which results 
in the power-disorder shock phase diagram of Fig. [T2{h). 
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FIG. 12: (Color online) (a)-(f) are the same as Fig. [TT] but 
at higher dye concentration (crhb = O.lmM); (g) the angular 
aperture 9 vs P: at this higher dye concentration the curves 
have recovered the threshold behavior also in the presence of 
disorder; (h) power-disorder phase diagram from the curves 
of(g)- 



VIII. CONCLUSION 

We have reported a detailed analysis of our experi- 
ments aimed at understanding the role of disorder in the 
occurrence of dispersive shock waves in a thermal defo- 
cusing medium. We collected the propagating and the 
transmitted intensity profile of a CW laser beam imping- 
ing on aqueous solutions of Rhodamine B, with an added 
controllable amount of disorder achieved by dispersing 



silica beads at well-defined concentrations. Resorting to 
the hydrodynamical approximation we analyzed the col- 
lected intensity distributions associated to the two ob- 
servables of the system: the shock point from the propa- 
gating intensity profiles and the angular aperture of the 
transmitted intensity profiles. Both the observables have 
evidenced the expected thresholds for the occurrence of 
the shock phenomenon with respect to the degree of non- 
linearity and in the amount of disorder. The calculation 
of the shock point has in fact led to the first determi- 
nation of disorder-power shock phase-diagram; also the 
trend of the angular aperture versus the laser power for 
the different silica concentrations has allowed to the cal- 
culation of two distinct shock diagrams related to the 
ordered and disordered cases. We also analyzed the de- 
gree of correlation of the shock images when increasing 
disorder. These experiments open the way to further 
investigations concerning the interplay between disorder 
and nonlinearity, with ramifications in several research 
directions, from basics physics, as the study of nonlin- 
ear waves in random media, to applied research, where 
the exploitation of nonlinear effects in disordered media, 
such as biological tissue and atmosphere, should be fun- 
damental in order to improve spectroscopy and imaging. 
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